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Abstract. Emission from the 3 Pi — 3 Po [Cl] transition 
at 492 GHz has been detected towards the dark cloud 
D478 in the Local Group galaxy M 31. Using literature 
detections of the lower 12 CO and 13 CO transitions, mod¬ 
els for the gas distribution in D 478 are discussed. The 
observed CO and C line ratios can be explained by two- 
component models (high density cores and low-density 
envelopes); single-density models appear less likely. The 
models indicate kinetic temperatures of the order of Tk « 
10 K. The beam-averaged column density of neutral car¬ 
bon is 0.3-0.8 times that of CO, whereas the total carbon 
to hydrogen ratio Nc/Nh is 5-3 x 10 -4 . The resulting CO 
to H 2 conversion factor X is about half that of the So¬ 
lar Neighbourhood. With temperatures of about 10 K and 
projected mass-densities of 5-10 M 0 pc -2 there appears 
to be no need to invoke the presence of very cold and very 
massive clouds. Rather, D 478 appears to be comparable 
to Milky Way dark cloud complexes at the higher rnetal- 
licity expected from its central location in M 31. In partic¬ 
ular, several similarities between D 478 and the Galactic 
Taurus-Auriga dark cloud complex may be noted. 
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1. Introduction 

Several molecular cloud complexes have been detected in 
the spiral arms of the Local Group galaxy M 31 by their 
12 CO <7=1-0 emission (Blitz 1985; Ryden & Stark 1986; 
Vogel et al. 1987; Casoli et al. 1987; Nakano et al. 1987; 
Lada et al. 1988; Casoli & Combes 1988; Berkhuijsen et 
al. 1993; Wilson & Rudolph 1993; Loinard et al. 1996b). 
Most appear to be giant molecular cloud complexes simi¬ 
lar to those in the disk of the Milky Way (Vogel et al 1987; 
Lada et al. 1988; Wilson & Rudolph 1993). However, Blitz 
(1985) noted velocity widths systematically higher than 
expected on the basis of the observed integrated line in¬ 
tensities. He proposed that the observed CO emission is 
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from both GMC’s and small clouds, and dominated by 
the contribution from the latter. Maps of a northeastern 
spiral arm of M 31 by Casoli et al. (1987) and Casoli & 
Combes (1988) and a southwestern spiral arm by Kutner 
et al. (1990) provide evidence for a dual population: in ad¬ 
dition to major complexes that have the characteristics of 
‘normal’ GMC’s, additional CO emission appears to orig¬ 
inate from numerous small clouds. Unbiased CO surveys 
along the minor and major axes of M 31 have been carried 
out by Sandqvist et al. (1989) and Loinard et al. (1995) 
respectively. The former found relatively strong CO emis¬ 
sion associated with M 31 spiral arm ‘4’ and weaker CO 
emission from the outer warp. The latter found CO emis¬ 
sion in the majority of positions sampled, in a number of 
cases without a counterpart dust cloud. They also found 
generally low 12 CO J=2-l/J=l-0 ratios which they in¬ 
terpreted as arising from cold (Tk < 5 K) molecular clouds 
at low density (nH 2 ~ 100 cm -3 ). 

In addition, they observed several dust clouds (Hodge 
1980) in the inner parts of M 31 (Allen & Lequeux 1993; 
Loinard et al. 1996a). Two of these, D 268 and D 478 were 
observed in the lower two transitions of both 12 CO and 
13 CO and subsequently modelled in detail by Allen et al. 
(1995). They conclude that the 12 CO emission is dom¬ 
inated by a very cold, low-density gas, while the 13 CO 
emission comes largely from higher-density clumps inside 
the clouds. In an attempt to further investigate the un¬ 
usual physical conditions of these molecular clouds, we 
have observed the strongest of these two objects, D 478, 
in the neutral carbon line. While this work was in progress, 
an interferometric map of the 12 CO <7=1-0 distribution, 
and the 12 CO 7=3-2 profile of D 478 have also become 
available (Loinard & Allen 1998). 

2. Observations 

The observations were carried out with the 15m James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea (Hawaii) 
[]. The 3 P! — 3 P 0 [Cl] transition at v = 492.161 GHz was 
observed towards D 478 for a total integration time of 100 

1 The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Particle Physics 
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Fig. 1. [Cl] spectrum observed towards D 478 in M 31. 
The vertical scale is = 0.53 T m b in K; the horizontal 
scale is velocity Vlsr in kms -1 . 

Table 1. [Cl] and CO in D 478 


Transition 

Tmb 

AV 

/ TnbdU 

Uhel 


(mK) 

kms -1 

K kms -1 

kms -1 

3 Pi - 3 P 0 [Cl] 

85 ±12 

15 ±2.5 

1.32 ±0.18 

-85.1 

12 CO ,7=2-1“ 

310 

19.5 

6.44 ±0.15 

-85.2 


Notes: a. Values in a 13” beam (J. Lequeux, private communi¬ 
cation) 

minutes (on+off) on November 29, 1996. At the observing 
frequency the resolution was 10.2" (HPBW) compared to 
a pointing accuracy better than 2" r.m.s. Weather condi¬ 
tions were excellent, resulting in an overall system tem¬ 
perature including sky of T sys = 1365 K. For a backend, 
we used the DAS digital autocorrelator system in a band 
of 250 MHz. The resulting spectrum was binned to a ve¬ 
locity resolution of 4.6 kms -1 . At this velocity resolution, 
the r.m.s. noise is about 14 mK in T£. We applied a lin¬ 
ear baseline correction only and scaled the spectrum to a 
main-beam brightness temperature, T m b = T^/r / m b using 
Tjmb = 0.53. Line parameters were determined by gaussian 
fitting; the results are given in Table 1. The spectrum is 
shown in Fig. 1. 

3. Results and analysis 


and Astronomy Research Council of the United Kingdom, the 
Netherlands Organisation for Scientific Research, and the Na¬ 
tional Research Council of Canada. 


3.1. [CIJ/CO ratio 

The 10" [Cl] profile in Fig. 1 is very similar to those of 
12 CO <7=2-1 and J=3-2 at 15" (Loinard & Allen 1998), 
suggesting similar spatial distributions of CO and C. The 
effect of resolution may be estimated by comparing the 
12 CO ,7=2-1 at a full resolution of 13" (kindly supplied 
by J. Lequeux, see Table 1) with the 12 CO result con¬ 
volved to 23" given by Allen et al. (1995). In amplitude, 
the full resolution profile is a factor 1.35 higher and in area 
it is a factor of 1.18 greater. We therefore feel confident 
in comparing the [CI]/ 12 CO <7=2-1 ratios resulting from 
our Table 1 to the CO ratios implicit in Allen et al. (1995 
- their Table 1). It then follows that the [Cl] 492 GHz 
line is similar in strength to the 13 CO <7=1-0 line, and 
considerably stronger than the 13 CO <7=2 1 line. Specifi¬ 
cally, we find [Cl]/ 13 CO(2-l) = 2.0 ±0.5. This is similar 
to the ratios 0.8 - 1.7 found at various translucent regions 
in the dark Galactic cloud L 183 (Stark et al. 1996), some¬ 
what higher than the corresponding range 0.9 - 1.3 found 
in the Galactic dark cloud TMC-1 (Schilke et al. 1995) 
and rather higher than the ratio of 0.8 characterizing the 
Orion Bar region (Tauber et al. 1995). For this reason, we 
decided to also include translucent conditions, character¬ 
ized by low CO and C column densities, in our analysis in 
addition to two-component models of the type considered 
by Allen et al. (1995), or the extremely low temperature 
model by Loinard & Allen (1998). 

3.2. Model calculations 

Our approach is the following. We first use radiative trans¬ 
fer models to explore the parameter space allowed by the 
observed line ratios. We then use the [Cl] observations and 
chemical model results (van Dishoeck & Black 1988; see 
also van Dishoeck 1998) to further constrain the parame¬ 
ters thus found. 

3.2.1. Radiative transfer modelling 

In order to constrain the range of possible cloud parame¬ 
ters consistent with the observed line ratios, we have per¬ 
formed statistical equilibrium calculations to determine 
the population distribution over the ground triplet levels 
of [Cl] and the lower 12 CO and 13 CO rotational levels, us¬ 
ing an escape probability method for the radiative transfer 
(cf. Jansen 1995; Jansen et al. 1994). These calculations 
include collisional excitation and de-excitation as well as 
spontaneous and stimulated radiative transfer. The equa¬ 
tions have been solved to find the level populations for 
a range of carbon and carbon monoxide column densi¬ 
ties N( 12 CO), N( 13 CO) and N(C), kinetic temperatures 
Tkin and H 2 volume densities n(H 2 ). We included a cos¬ 
mic background radiation field of 2.75 K. and an incident 
radiation held 7uv ~ 0.5 (7uv corresponding to 7iooo = 
4.5 x 10 -8 photons s -1 cm -2 ). 
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Table 2. Model parameters and line ratios in D 478 




Single-Component Models 

Two-Component Models 

Observed 

Tkin 

(K) 

6 

10 

14 

6, 6 

6, 10 

10/10 

10/14 


Weight 


— 

— 

— 

1/2, 1 

1, 1 

1/6, 1 

1/4, 1 


n(H a ) 

(cm -3 ) 

5000 

1500 

750 


3000, 

100 



N( 13 CO)/dV (10 14 cm" 

'7 kms -1 ) 

4 

7 

8 

20, 2 

20, 2 

60, 6 

52, 5 


N( 12 CO)/dF (10 16 cm" 

'7 kms -1 ) 

0.8 

1.3 

1.7 


72, 7.2 



N{C)/dV (10 16 cm - 

' 2 / kms -1 ) 

6 

2.5 

2 

120, 160 

30, 105 

15, 80 

5, 70 


12 CO (2—1)/(1—0) 


0.46 

0.46 

0.46 

0.40 

0.46 

0.45 

0.48 

0.42 ±0.08“ 

12 CO (3-2)/(2—1) 


0.12 

0.21 

0.26 

0.32 

0.30 

0.28 

0.32 

0.33 ± 0.10 6 

13 CO (2—1)/(1—0) 


0.40 

0.40 

0.40 

0.34 

0.35 

0.48 

0.50 

0.45 ±0.15“ 

12 CO/ 13 CO (1-0) 


7.4 

7.8 

7.8 

7.4 

7.4 

8.5 

8.7 

8.8 ± 1.8“ 

12 CO/ 13 CO (2-1) 


8.5 

8.7 

9.1 

8.4 

9.7 

7.5 

8.5 

8.4 ± 1.7“ 

CI/ 12 CO (2-1) 





0.24 




0.24 ± 0.05“ 


Notes: a. Derived from Allen et al. (1995); b. Derived from Loinard & Allen (1998); c. This Paper. 


As we have four independent CO line ratios, single- 
component CO models are, in principle, fully constrained 
by the four adjustable parameters. In practice, relatively 
large observational uncertainties allow only fitting of three 
parameters as a function of the fourth, for which we have 
chosen Tki n = 6, 10 and 14 K. Neutral carbon column den¬ 
sities are adjusted to yield, at the relevant kinetic temper¬ 
ature, the observed [Cl]/ 12 CO (2-1) ratio of 0.24. By com¬ 
bining these with the model calculations of van Dishoeck 
& Black (1988), we then estimate the gas-phase carbon 
fraction N C /N H = (N(CO) + N(C))/(2N(H 2 + N(HI)). 
From Bajaja & Shane (1982) and Brinks & Shane (1984) 
we find a neutral hydrogen column density N(HI) = 1.8 
x 10 20 cm -2 at the position of D 478. 

Parameters thus determined for single-component 
models are shown in Table 2, together with corresponding 
model line ratios. Although these models generally yield 
greater differences between the 12 CO (2-1)/(1-0) and (3- 
2)/(2-1) ratios than observed, the 14 K model fits the 
observations within the errors. The 10 K model is barely 
consistent and the 6 K model fails to reproduce the ob¬ 
served ratios, predicting too low a 12 CO J = 3-2 strength. 
Loinard & Allen (1998) have argued that the 12 CO struc¬ 
ture of D 478 can be explained by a large, dense and very 
smooth cloud in LTE at even lower temperatures close 
to that of the cosmic background. Although it is possi¬ 
ble to reproduce the 12 CO line ratios in this manner, the 
observed 13 CO (2—1)/(1—0) ratio of about 0.5 is only ob¬ 
tained at high 13 CO optical depths requiring 12 CO/ 13 CO 
ratios of order unity, which is clearly not the case. Alterna¬ 
tively, the observed 12 CO/ 13 CO ratios can be reproduced 
with lower 13 CO optical depths satisfying one of the two 
observed 12 CO/ 13 CO ratios, but the 13 CO (2—1)/(1—0) 
ratio then is only a third of the actually observed value. 

For such reasons, Allen et al. (1995) concluded to the 
necessity of a model incorporating at least two compo¬ 
nents. In their analysis, they modelled the CO emission 


from D 478 with two different density/temperature com¬ 
ponents, and showed that the observed line ratios are, for 
instance, reproduced by a cold low-density gas (Tk « = 4 
K, = 100 cm -3 , dV = 5 kms -1 ) and a high-density 
gas (Tk = 7 K, Tin = 3000 cm -3 , dV = 10 kms -1 ), both 
with a low surface-filling factor. 

We have also calculated two-component models. As 
there are eleven adjustable parameters (six column densi¬ 
ties, two densities, two temperatures and the relative con¬ 
tributions of the two components to the total emission), 
no unique solution can be obtained. We have chosen to 
model the same H 2 densities as given in the example by 
Allen et al. (1995) and similar 12 CO column densities, 
and calculate the remaining parameters as a function of 
temperature. Furthermore, we require the column density 
ratio of the tenuous and the dense components to be the 
same for 13 CO and for 12 CO. Over the range of appli¬ 
cable parameters, the [Cl] emission from the dense com¬ 
ponent is more efficiently produced (i.e. characterized by 
a higher ratio of [Cl] intensity to carbon column density) 
than that from the tenuous component. By taking the 
maximum permissible contribution from the dense com¬ 
ponent, and ascribing the remainder to the tenuous com¬ 
ponent, we have minimized overall neutral carbon column 
densities required, and thus maximized permitted total 
hydrogen column densities. Although other combinations 
of input parameters are possible, the range of admissible 
parameters is nevertheless limited by the observations. For 
instance, changes in the assumed H 2 densities and 12 CO 
column densities, immediately require adjustment of the 
relative weights of the two components in order to repro¬ 
duce the same line ratios which tends to minimize the 
effect of these changes on the parameter ratios as well. 

In Table 2, we present four such cases, for tempera¬ 
tures varying from 6 to 14 K. As an example, the 6/6 K 
case corresponds closest to the example given by Allen et 
al. (1995), with the surface filling factor of the low-density 
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component twice that of the high-density component. For 
this case, we obtain a reasonable fit for a 12 CO/ 13 CO 
column-density ratio of 360. If, instead, we require the 
column density ratio to be 90, 13 CO optical depths in¬ 
crease significantly and we obtain for the same relative 
filling factors a 12 CO/ 13 CO temperature ratio of less than 
4, inconsistent with the observed values. This may be re¬ 
paired by increasing the relative filling factors to 1:6 or 
higher, but then the 12 CO J=2-l/J=l-0 and especially 
the J=3-2/J=2-l ratios drop to unacceptably low values 
of 0.32 and 0.19 respectively. 

As Table 2 shows, increasing the temperature generally 
leads to higher 12 CO/ 13 CO ratios, lower ratios of neutral 
carbon to CO column density and a lower required contri¬ 
bution of the dense component. The dense component has 
low N(C)/N(CO) ratios characteristic of Giant Molecu¬ 
lar Cloud complexes in the Milky Way (0.1-0.2: Keene et 
al. 1985), whereas the tenuous component exhibits ratios 
similar to those found in Galactic translucent and dark 
clouds (0.3-3: Stark & van Dishoeck 1994; Stark et al. 
1996). 

If we increase the temperature of the dense compo¬ 
nent over that of the tenuous component, the 13 CO col¬ 
umn densities change only by a small amount, the implied 
C column densities decrease and the relative contribution 
of the tenuous component also decreases. The predicted 
line ratios, although somewhat different, are consistent 
with the observed ratios. Thus, from a radiative-transfer 
point of view, the temperature of D 478 is not strongly 
constrained in a two-component model, although we note 
that the higher temperature (10/10 K and 10/14 K) mod¬ 
els give slightly better fits to the observed line ratios. In 
the models in Table 2, the [Cl] emission is optically thick 
at 6 K, and optically thin at 14 K. In all cases, the ex¬ 
pected strength of the 809 GHz line is (much) less than 
10 per cent of that of the 492 GHz line. As neutral car¬ 
bon radiates ineffectively at temperatures below 10 K, the 
observed [Cl] intensity rapidly requires very large C col¬ 
umn densities and high N(C)/N(CO ) ratios if lower tem¬ 
peratures are assumed. This again rules out satisfactory 
solutions at very low temperatures. 

4. Discussion 

4-1 ■ The environment of D 478 

We have used the chemical models by van Dishoeck & 
Black (1988) assuming an incident radiation field 7uv ~ 
0.5, corresponding to Ziooo = 2.25 x 10 -8 photons s _1 
cm -2 . Is this a realistic value? No measurements at A = 
100 nm exist for M 31, but we may estimate the value of 
7uv at the location D 478 by using data obtained at 155 
nm (Carruthers et al. 1978; Wu et al. 1980; Israel et al. 
1986). After a correction of E(B — V) = 0.11 for Galactic 
foreground extinction only, we find from these data an ir¬ 
radiation of D 478 by about /1550 = 2.5 x 10 -8 photons 


s ^ 1 cm -2 . A similar result was obtained independently by 
Koper (1993). About 60% of this originates in the rela¬ 
tively nearby M 31 bulge, while about 40% is contributed 
by the more distant star-forming ring. A somewhat un¬ 
certain extrapolation of the 155 nm fluxes to A = 100 nm 
then suggests irradiation by /yv ~ 0.3-0.4, i.e close to 
the assumed value. A substantially lower value of /yv ~ 
0.1 would change our conclusions only slightly: it would 
require lower H 2 column densities, and somewhat higher 
carbon abundances. In the absence of measurements, we 
have neglected the possible presence of ionized carbon 
(C+). If C + is present in any significant amount, it would 
raise the total carbon column Nq. and lower the fraction 
of all carbon contained in CO. In turn, this would tend 
to lower the resulting H 2 column densities, and lead to 
significantly higher carbon abundances. 

The models would be better constrained if we had an 
independent way of estimating total hydrogen or molec¬ 
ular hydrogen column densities. This would be the case 
if, for instance, the CO intensity to H 2 column density 
conversion factor X or the carbon abundance in the cen¬ 
ter of M 31 were known. Unfortunately, neither is well- 
determined. In the disk of our Galaxy, a conversion fac¬ 
tor of X = 2 x 10 2 ° mol cm -2 (Kkms -1 ) -1 (Strong & 
Mattox 1996) applies, within a factor of two. However, 
in the central part of the Galaxy X is a factor of 3-10 
lower (Sodroski et al. 1995), whereas in low-metallicity 
environments X may be up by one or even two orders of 
magnitude (Israel 1997). 

Oxygen abundances in HII regions and supernova rem¬ 
nants in the disk of M 31, studied by Blair et al. (1982), 
indicate the presence of a gradient, whose extrapolation 
inwards yields an [O]/[H] abundance of about 1.3 x 10 -3 
at the galactocentric radius of D 478. Data summarized by 
Garnett et al. (1995) and Kobulnicky & Skillman (1998) 
suggest that at such high oxygen abundances, the ratio 
of carbon to oxygen abundances is given by log [C]/[0] 
= -0.2 ± 0.3 so that [C]/[H] = 8(+8, -4) x 1CT 4 . This 
result has a large uncertainty, however, because no direct 
abundance measurements exist within R = 5 kpc, because 
the M 31 disk data exhibit a relatively large scatter, and 
because the relation between [0]/[H] and [C]/[0] abun¬ 
dances is not well-established at such high [0]/[H] abun¬ 
dances. In addition, a significant fraction of all carbon in 
quiescent dark clouds such as D 478 is expected to be de¬ 
pleted onto dust grains, leaving actual gas-phase carbon 
abundances substantially lower, by factors of order <5c ~ 
0.4. This would lead us to expect total column density 
ratios Nq/Nh «3x 10~ 4 . 

Moreover, an inward extrapolation from the disk may 
not be representative of conditions in D 478. The velocity 
of D 478 is anomalously low for its projected distance to 
the nucleus. In position-velocity maps presented by Brinks 
& Shane (1984), one of which is reproduced by Loinard et 
al. (1995), it coincides with a loop-like HI structure. This 
structure exhibits a north-south symmetry with respect 
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Table 3. Beam-averaged parameters for D 478 


Units Single-Component Models Two-Component Models 


Tkin 

(K) 

6 

10 

14 

6, 6 

6, 10 

10, 10 

10, 14 

N( 12 CO) 

(10 16 cm” 2 ) 

5 

4 

4 

160 

135 

50 

45 

N( C) 

(10 16 cm” 2 ) 

26 

6 

4 

130 

35 

25 

15 

A(C)/N(CO) 


5 

1.5 

1 

0.8 

0.3 

0.5 

0.3 

N c 

(10 16 cm” 2 ) 

31 

10 

8 

290 

170 

75 

60 

Nn 

(10 21 cm” 2 ) 

1.6 

1.7 

1.2 

3.5 

4 

2 

2 

N c /N h 

(10- 4 ) 

2 

0.6 

0.7 

8 

4.5 

3.5 

3 

N{ H a ) 

(10 21 cm” 2 ) 

0.7 

0.7 

0.5 

1.7 

1.8 

1.0 

0.9 

/(H 2 ) (dense) 


0.22 

0.11 

0.10 

0.27 

0.27 

0.06 

0.06 

/(H 2 ) (tenuous) 


— 

— 

— 

0.54 

0.27 

0.33 

0.24 

A( 12 CO)/A(H 2 ) 

(10- 4 ) 

0.7 

0.5 

0.8 

9.5 

7.5 

5 

5 

dM/dA 

(M 0 pc 2 ) 

4 

4 

3 

9 

9 

5 

5 

X 

(10 20 cm” 2 (Kkins” 1 )” 1 ) 

0.5 

0.5 

0.3 

1.1 

1.2 

0.7 

0.6 


to the center, but is most pronounced north of the cen¬ 
ter. It appears to represent rapidly rotating gas in the 
inner part of M 31, moving in inclined elliptical orbits in 
a weak bar-like potential (cf. Stark & Binney 1994). D 478 
is located close to the northern tangential point of this ro¬ 
tating structure. It is conceivable that this is the signature 
of a merging event in which a small late-type galaxy fell 
into the large spiral M 31. Such a notion is supported by 
the presence of a double nucleus in M 31 (cf. Lauer et al. 
1993). In that case, the metallicity of the inner gas clouds 
may be substantially lower than expected from extrapola¬ 
tion of the disk gradient. If the metallicity is comparable 
to that of the outer disk of M 31 we have [O] / [H] = 4 x 
10” 4 and again following Garnett et al. (1995) we find a 
low total carbon abundance [C]/[H] = 1.5 x 10” 4 . Taking 
into account depletion, this case would lead us to expect 
N c /N n « 0.7 x 10” 4 . 

4-2. Conditions in D 478 

In Table 3, we have used the ratio of the observed to 
model line intensities to calculate beam-averaged column 
densities of carbon monoxide and neutral carbon, as well 
as those of total carbon (C + CO), total hydrogen, and 
molecular hydrogen. We have used the total hydrogen col¬ 
umn densities, multiplied by 1.35 to take into account a 
contribution by helium, to calculate the projected mass 
density dM/dA (corrected by cos?' with i = 77°) and we 
also determined the CO - H 2 conversion factor X. As 
shown in the preceding section, we have, if anything, over¬ 
estimated hydrogen column densities (and masses), and 
underestimated carbon abundances. In addition, we have 
estimated surface filling factors /(H 2 ) by dividing the ob¬ 
served velocity-integrated line intensities by the model in¬ 
tensities, and assuming that the surfaces covered by CO 
and by H 2 are identical. 


First we consider the single-density models. The phys¬ 
ical equivalent of these models is that of high-density fil¬ 
amentary clouds rather uniformly but incompletely fill¬ 
ing a large surface area. We have already ruled out the 
model for T k = 6 K, because it does not reproduce the 
observed 12 CO J= 3-2 intensity. Both the 10 K and 14 K 
models fit the observed intensities although the predicted 
12 CO J= 3-2 intensity is low compared to the observed 
value (Table 2), rendering the 10 K model in particular 
somewhat marginal. However, single-component models 
are consistent only with low [C] / [H] ratios. They are ruled 
out if we accept the extrapolated high carbon abundances, 
but would be just consistent with an anomalously low 
metallicity as discussed above. The single-density mod¬ 
els furthermore require rather low 12 CO/ 13 CO abundance 
ratios of about 20. This is not out of the question, because 
in cold diffuse and translucent clouds, such as described 
by the single-density models in Table 2, ion-molecule ex¬ 
change reactions may favour the formation of 13 CO at 
the expense of 12 CO for temperatures T k < 36K although 
13 CO remains more susceptible to photodissociation than 
12 CO. Surface filling factors are between 0.2 and 0.1, the 
lower values pertaining to the higher temperatures. The 
relatively low H 2 column densities imply projected mass 
densities of only about 4 M 0 pc -2 (Table 3), and X values 
about five times lower than those in the Solar Neighbour¬ 
hood. We conclude that for D 478 single-component mod¬ 
els at temperatures T k i n > 10 K on the one hand cannot 
be ruled out, but on the other hand do not appear to be 
very convincing. 

The two-density models provide better fits to the ob¬ 
servations; like the single-component models they have av¬ 
erage H 2 column densities below the star formation tresh- 
old in keeping with the quiescent appearance of D 478. 
The 6/6 K model requires a relatively large Nc/N-n ra¬ 
tio of eight, implying a total carbon abundance [C]/[H] 
> 2 x 10” 3 which seems rather high. More reasonable 
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carbon abundances would imply higher H 2 column densi¬ 
ties and projected gas mass densities, but the large neu¬ 
tral carbon column densities required by the observed line 
strength at the assumed temperature of 6 K would then 
imply N(C)/N(CO) ratios much larger than allowed at 
such column densities. 

With increasing temperature, the CO and [Cl] column 
densities required by the observed line emission become 
more modest. Because the carbon/hydrogen column den¬ 
sity ratio also drops, total hydrogen and molecular hy¬ 
drogen column densities change relatively little. At these 
temperatures, Nq/Nh «4x 10” 4 , which corresponds to 
a total carbon abundance [C]/[H] «lx 10” 3 if 5c = 0.4. 
This is about 2.5 times the value ascribed to the Solar 
Neighbourhood, and close to the value expected from the 
extrapolation of M 31 disk abundances. The CO inten¬ 
sity to H 2 column density value for these models is rela¬ 
tively well-determined with X = 0.9±0.3 xlO 20 mol cm” 2 
(Kkms” 1 )” 1 , or about half of the Solar Neighbourhood 
value. Projected mass densities are in the range of 5-9 M 0 
pc” 2 . Surface filling factors for the emission within the 
beam are given in Table 3 for both the high-density and 
the low-density components; the ratio of the two is for each 
model, of course, equal to the ratio of the weights given 
in Table 2. Again, the lower temperature models have the 
largest filling factors. The two models with dense gas tem¬ 
peratures of 10 K yield surface filling factors of 6% for the 
dense gas and about 30% for the tenuous gas. These val¬ 
ues appear to be consistent with the small fraction (15%) 
of the single-dish CO flux recovered in the interferometric 
map by Loinard & Allen (1998). The model H 2 column 
densities divided by the input H 2 densities provides a mea¬ 
sure for the line-of-sight depth z through the clouds in the 
emitting surface fraction. For the two-component models 
in Table 3, we find 0.3 pc j Zd j 0.7 pc and 4.5 pc j Zt j 
7.5 pc. The subscripts d and t refer to the dense and the 
tenuous components respectively, and the larger depths 
occur at the higher temperatures. Both the surface filling 
factors and the subcloud depths suggest that the structure 
of D 478 is highly fragmented and filamentary. 

The models discussed do not fully constrain the phys¬ 
ical conditions applying to D 478. Various other combina¬ 
tions of parameters could be used as input. Multi-density 
models incorporating temperature gradients might be bet¬ 
ter representations of the actual structure of D 478, and 
could easily be made to satisfy the observations. For in¬ 
stance, adding a contribution of material at intermediate 
densities and temperatures to the two-component mod¬ 
els in Table 2 within the limits imposed by observed line 
ratios would indeed change the detailed set of physical pa¬ 
rameters, but would not substantially change the overall 
result. We conclude that the observations of D 478, and 
in particular the relatively strong [Cl] emission appear to 
rule out the presence of very large amounts of molecular 
gas at very low temperatures in this complex. At the very 
least, our models show that there is no compelling need to 


assume such a state of affairs. The interferometric results 
obtained by Loinard & Allen (1998) indicate a relatively 
small contribution by dense, compact material. This sug¬ 
gests that the 10/10 K and 10/14 K models with their 
small surface-filling fraction of the dense component are a 
better representation than the lower-temperature models. 

Most likely, D 478 is rather similar to dark cloud com¬ 
plexes in our own Galaxy, at the higher metallicity con¬ 
sistent with its central location in M 31. It has a modest 
projected mass density and a relatively low V-factor, both 
of which are consistent with this conclusion. A comparison 
with the Galactic Taurus-Auriga dark cloud complex at a 
distance of 140 pc is illustrative. We have already noted 
similar [CI]/ 13 CO (2-1) ratios (Schilke et al. 1995). Like¬ 
wise, the Taurus 12 CO (2 1)/(1-0) ratio of 0.53 (Sakamoto 
et al. 1997) is close to that of D 478. More in detail, an area 
30° (80 pc) across was mapped in CO by Ungerechts & 
Thaddeus (1987). Ignoring the more distant and unrelated 
clouds, we find that the mean integrated line intensity of 
the complex is 4.5 Kkms” 1 with individual clouds having 
intensities between 3 and 15 Kkms” 1 . The surface filling 
factor of the complex is 0.3, and its appearance is indeed 
very fragmented and clumpy. The complex has a mean 
mass-density dM/dA = 7 M Q pc” 2 . Both the observed 
line ratios and the derived quantities thus suggest that 
D 478 is similar to a larger version of the Taurus-Auriga 
dark cloud complex. 

Assuming an overall cloud size of 1' (Allen & Lequeux 
1993) and including a contribution by helium, the gas 
mass of D 478 is of the order of 1 x 10 6 Mq. This is an 
order of magnitude lower than the virial mass deduced by 
Allen & Lequeux (1993) and Loinard & Allen (1998). How¬ 
ever, it is doubtful whether molecular cloud complexes are 
virialized at all. Moreover, the shape of the line profiles 
of D 478 suggests the presence of at least two distinct 
components. Its location near the tangential point of a 
rapidly rotating structure with noncircular motions sug¬ 
gests that various individual clouds, separated by signifi¬ 
cant distances, may be seen with very similar velocities in 
the same line of sight; the appearance of the HI distribu¬ 
tion in the position-velocity maps published by Brinks & 
Shane (1984) lends credibility to this suggestion. In this 
respect it is relevant to note that the velocity width of line 
spectra incorporating emission of various unrelated clouds 
in the same beam will reflect the cloud-cloud velocity dis¬ 
persion of about 9 kms” 1 (Stark 1979) in addition to any 
systematic velocity shifts and the velocity widths of the 
individual clouds. Applying the virial equation to such a 
collection of clouds will lead to a potentially large overes¬ 
timate of the actual mass. 

5. Conclusions 

1. We have detected the 492 GHz emission line from the 

3 Pi — 3 Po [Cl] transition towards the dark cloud D 478 

in M 31. The [Cl] line strength is comparable to that of 
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13 CO <7=1—0 and about twice as strong as that of the 
13 CO ,7=2-1 transition. This is similar to the intensity 
ratios found in Galactic dark and translucent clouds. 

2. We have used 12 CO and 13 CO intensities from the lit¬ 
erature together with the newly determined [Cl] inten¬ 
sity to model conditions in the D 478 cloud complex. 
We have considered both single-density filamentary 
clouds and two-density core/envelope clouds. In partic¬ 
ular, inhomogeneous (two-component) models produce 
good fits to the observations. Dense core N ( C)/N(CO) 
ratios are found to be similar to those of Giant Molec¬ 
ular Clouds in the Milky Way, whereas the ratios de¬ 
rived for the extended more tenous gas are similar to 
those found in Galactic dark and translucent clouds. 

3. D 478 appears to be similar to Galactic dark cloud 
complexes, such as the Taurus-Auriga dark cloud com¬ 
plex. Most likely, it is characterized by a higher metal- 
licity consistent with an extrapolation of M 31 disk 
abundances to its central location in the galaxy. Its 
kinetic temperature is of the order of 10 K. There ap¬ 
pears to be no need to assign very low temperatures 
and very high H 2 mass densities to the D 478 complex; 
our results do not support the virial mass surface den¬ 
sities of D 478 of ~100 M(7) pc -2 suggested by Loinard 
& Allen (1998). 

4. Actual cloud mass densities projected onto the M 31 
midplane are of the order of 5-10 Mgpc -2 . The CO 
to H 2 conversion factor is X = 0.9 ±0.3 x 10 2 ° cm -2 
(K km s~ 1 ) “ 1 , about half the value in the Solar Neigh¬ 
bourhood. 

5. The D 478 complex may not be a single entity, but 
instead may consist of various clouds at different lo¬ 
cations projected onto the same line-of-sight towards 
the tangential direction of a rapidly rotating structure 
characterized by noncircular velocities. 
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